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Abstract 

Averaging is an important method to extract effective macroscopic 
dynamics from complex systems with slow modes and fast modes. 
This article derives an averaged equation for a class of stochastic par- 
tial differential equations without any Lipschitz assumption on the 
slow modes. The rate of convergence in probability is obtained as a 
byproduct. Importantly, the deviation between the original equation 
and the averaged equation is also studied. A martingale approach 
proves that the deviation is described by a Gaussian process. This 
gives an approximation to errors of 0{e) instead of 0{^/e) attained in 
previous averaging. 

Keywords: Slow-fast stochastic partial differential equations, averaging, 
martingale 



1 Introduction 

The need to quantify uncertainties is widely recognized in modeling, analyz- 
ing, simulating and predicting complex phenomena [Sllinillll e.g.]. Stochas- 
tic partial differential equations (spdes) are appropriate mathematical mod- 
els for many multiscale systems with uncertain influences [17]. 
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Very often a complex system has two widely separated timescales. Then 
a simplified equation which governs the evolution of the system over the 
long time scale is highly desirable. Such a simplified equation, capturing the 
dynamics of the system at the slow time scale, is often called an averaged 
equation. There is a great deal work on averaging principles for determinis- 
tic ordinary differential equations PQ [21 [ISl e.g.] and for stochastic ordinary 
differential equations [HI El El e.g.]. But there are few results on the aver- 
aging principle for SPDEs. Recently, an averaged equation for a system of 
reaction-diffusion equations with stochastic fast component was obtained by 
a Lipschitz assumption on all nonlinear terms [1]. The resultant averaged 
equation is deterministic. 

This article derives an averaged equation for a class of SPDEs with stochas- 
tic fast component and proves a square-root rate of convergence in probabil- 
ity. Furthermore, the deviation between the original system and the averaged 
system is determined. 

Let D be an open bounded interval and L'^{D) be the Lebesgue space of 
square integrable real valued functions on D. Consider the following slow-fast 
system 

du' = [Au' + f{u\ v')] dt + (Ti dWi{t) , u'{0) = uq e L\D) , (1) 
dv' = -[Av' + g{u',v')]dt + ^dW2{t) , v' (0) = vq e L\D) , (2) 

with Dirichlet boundary condition. Here Wi and W2 are mutually inde- 
pendent L'^{D) valued Wiener processes defined on a complete probability 
space JF, P) detailed in the following section. If for any fixed -u, the fast 
system ([2]) has unique invariant measure /i^, then as e ^ , under some 
conditions, the solution of ([T]), converges in probability to the solution of 

du = [Au + f{u)]dt + aidWi{t) , (3) 

m(0) = Mo and u\dD = . (4) 

Here, the average 

fiu) = / f{u,v)nu{dv) . (5) 

JH 

And the convergence rate is proved to be 1/2 in the following sense for any 
K > 

P| sup \u'{t) - u{t)\ < C^t^'^] >l-K (6) 

for some positive constant C^; see Section HI 

Furthermore by estimate ([6]), as e — > the limit of {u'it) — u{t))/^ is 
proved to be a Gaussian process; see Section [5l 
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We stress that Theorem [S] gives a much better approximation than the 
averaged equation. 



2 Preliminaries and main results 

Let H = L'^{D) with L^-norm denoted by | ■ | and inner product by (■,■). 
Define the abstract operator A = A with zero Dirichlet boundary condition 
which defines a compact analytic semigroup e^*, t >0 on H. For any a > , 
u E H define \u\a = \A°'^'^u\ and for a = 1 , the norm is denoted as || ■ ||. Then 
let Hq be the space the closure of C^{D), the space of smooth functions with 
compact support on D, under the norm \ ■ \a ■ Furthermore, let H~°' denote 
the dual space of iig and denote by Ai the first eigenvalue of A . Also we 
are given H valued Wiener processes Wi{t) and H^2(^), ^ > 0, which are 
mutually independent on the complete probability space JF, JF^, P) [T2] . 
Denote by E the expectation operator with respect to P. Then consider the 
following SPDEs with separated time scale 



Here ai G M , a2 7^ are arbitrary real numbers. For our purpose we adopt 
the following four hypotheses. 

Hi f{x,y) : R X R — i> M is continuous, there is a positive constant Cf, such 
that f'^{x,y) < Cf , \ fy{x,y)\ < Cf and that for any x,y 



for some positive constants a, b and c. 

H2 g{x,y) : R X R — i> R is continuous and is Lipschitz with respect to the 
both variables with Lipschitz constant Cg. For any x,y 




(7) 
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\fix,y)\^<ax' + by^ + c, 

/(x, y)x < —ax^ — bxy + c , 

{f{xi,y) - f{x2,y)){xi - X2) < a{xi - X2f + c. 



9{x,y)y < -dy'^ + exy 



for some positive constants d and e. 



H3 b> e and < Ai . 
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H4 Wi and W2 are Q- Wiener processes with covariance operator Qi and Q2 
respectively. Moreover, tr[yl^/^(5j] < 00 , i = 1,2 . 

With the above assumptions we have the first result on the fast component 
that is proved at the end of Section El 

Theorem 1. Assume Hi H^. For any fixed u E H , system ^ has a unique 
stationary solution, rjl^if), with distribution /i^ independent of e. Moreover, 
the stationary measure is exponentially mixing. 

Then we prove the following averaging result. 

Theorem 2. Assume Hi H4. Given T > 0, for any uq E H , solu- 
tion u^it, Uq) of (Cy converges in probability to u in C(0, T; H) which solves 
Moreover the convergence rate is 1/2 that is for any k > 



I sup \u\t) - u{t)\ < C!^^/e] >1-K 



0<t<T 

for some positive constant > . 

Having the above averaging result we consider the deviation between u"^ 
and u. For this we introduce 

z^{t) = ^{u^-u). (9) 

Then we have 

Theorem 3. z"^ converges in distribution to z in space C{0,T]H) which 
solves 



z = Az + f'^{u)z+ ^B{u)W (10) 
where Biu] : H ^ H is Hilbert-Schmidt with 



Bin) = 2 E 



00 







ifiu,vu{t)) - /») {f{u,nM) - /») 



dt 



fM= / f:,iu,v)^,-{dv) 



H 



and W{t) is an H -valued cylindrical Wiener process with covariance opera- 
tor Idn- 
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3 Some a priori estimates 



This section gives some a priori estimates for the solution of (I71)-(IH1) which 
yields the tightness of u'' in space C(0,T; H). First we give a wellposedness 
result. 

Theorem 4. Assume Hi H4. For any uq (z H , vq & H and any T > 

0, there is unique solution {u'{t),v'{t)) in L'^{n,C{0,T; H) f] L'^{0,T; H^)) 
for 

The above result is derived by a standard method [12] and so is here 
omitted. Then we have the following estimates for {u^,v'^). 

Theorem 5. Assume Hi H4. For Uq G and Vq G , for any T > , 

there is a positive constant Ct which is independent of e, such that 

E sup ||M'(t)f + sup ¥.\v'{t)\^ <Ct(\\u4'^ + \vo\^) (11) 

0<t<T 0<t<T 

and for any positive integer m, 

e[ \\u'{s)f"'ds + E [ \\v'{t)fds<CT{\uo\^ + \vo\^). (12) 







Moreover 

nu'\mO,T;H-^) < CtHuoI^ + l^^oH • (13) 

Proof. Applying Ito formula to |M'^(t)p and |w'^(t)p respectively and by Gron- 
wall lemma there is positive constant C which is independent of e such that 
for any t > 

E|M^(t)P + eEK;^(t)p < C{\uo\^ + \vo\^) . (14) 

At the same time we have for any T > that there is positive constant Ct 
which is independent of e such that for < t < T 

e[ \\u'{s)\\^ds + E [ \\v'{s)\\^ds <Ct{\uo\^ + \vl\) (15) 
Jo Jo 

and 

supE|T;^(t)P < Ci\uo\^ + \vo\^). (16) 
t>o 

By applying Ito formula to ||'u'^(i(:)|p and using ( |T5|) . a lemma of Da Prato 
& Zabczyk's [T2| Lemma 7.2] proves 

E sup Wu^m^ <CT{\\uof + \vo\^). (17) 

0<t<T 

Further, applying Ito formula to ||M''(t)|p™' and ||f^(t)|p"*, by the assump- 
tion on / and ( JTSl) . we then establish the inequality (JT2l) . The bound ( JT3l) 
can be proved by estimating ( JT2l) and the embedding C L^{D). □ 
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Now we show that C^u"), the distribution of m", is tight in C{0,T; H). 
For this we need the following lemma by Simon 



Lemma 6. Let E, Eq and Ei be Banach spaces such that Ei (£ Eq , the 
interpolation space {Eq, Ei)e^i C E with 9 G (0, 1) and E C Eq with C 
and (£ denoting continuous and compact embedding respectively. Suppose 
Po,Pi G [1, oo] and T > , such that 

V is a bounded set in L^^(0,T; Ei) 

and 

dV := {dv.ve V} IS a bounded set m Lp°{0,T; Eq). 
Here d denotes the distributional derivative. If 1 — 9 > l/pg with 

1 1-9 9 



Pe Po Pi 

then V is relatively compact in C(0,T; E). 

By the above lemma we have the following result. 

Theorem 7. Assume Hi H4. {£(u'^)}e>o is tight in space C{0,T] H) . 

Proof. Taking Eq = H-\ E = H and Ei = Hi and po = 2 , ^ = 1/2. 
By Theorem [5], pi can be taken as arbitrary lager positive integer, then by 
Lemma El {£(M^)}e>o is tight in space C(0,T; H). □ 

Proof of TheoremUi For any two solutions vl and V2, the Ito formula yields 

E\vl{t) - < e-2(^i-^«)*/^E|t;i(0) - 



which means the existence of a unique stationary solution r^^ for ([H]) dis- 
tributes as fiu such that for any Vq E H 

nv'it) - < e-^(^-^«)*/^E|t;o - vim' , (18) 

which yields the exponential mixing. Moreover, since \fy{x,y)\ < Cf, we 
also have 

Ef{u,v'{t))- I fiu,x)fiu{dx) <C(l + |^;o^e-'(^^-^«)*/^ (19) 
Jh 

By the time scale transformation t ^ r = et , is transformed to 

dv = [Av + g{u, v)] dr + a2 dW2{T) , v{0) =vo, (20) 
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where W2 is the scaled version of W2 and with the same distribution. Then 
the SPDE (120]) has a unique stationary solution rju with distribution And 
by the ergodic property of we have 



1 f 

hm - / f{u,r]u{s)) ds . 



t^oo t 



(21) 







Furthermore, by a generalized theorem on contracting maps depending on a 
parameter [3l Appendix C], [3], f„(t) is differential with respect u with 



sup 

u,vo(iH, 0<t<oo 



\DuVu\c{H) < C 



for some positive constant C . 



[22) 



□ 



We end this section by giving the following a priori estimates on the 
solutions of the averaged equation ([3])-(|l]) . First we need a local Lipschitz 
property of / which is yielded by (1^ and the following estimate, for any 

ui , ui G Hi 



T 
1 

< - 

r 



< 2 



f «i (s)) - /(m2, f „2 (s))] ds 

[f{Ui, f (s)) - f{u2, Vu^ (s))] rfs 



H 



1 

r 

ll^iil 



[/(^^2,tV(-s)) - f{u2,Vu^{s))\ds 



H 



|m2||^+ sup \DuV{t)\c(H) 

u^VQdH, 0<t<oo 



Ml — U2\ ■ 



(23) 



Lemma 8. Assume Hi H:^ . For any Uq G H^, for any T > 0, ^-(^ 
has a unique solution u G L^(r2, C(0, T; H) fl L^(0, T; ifg)) . Moreover there 
is a positive constant Ct such that for any positive integer m > 2 and any 
0<t<T 

E\\uit)\r<CTil + \\uo\n- (24) 
Proof. Applying Ito formula to |M(t)p yields 

~\u{t)\' = -MtW + {f{u{t))Mt)) + Mu{t), Wi) + ^ trQi . 



By (EH) and assumption Hi 



{f{u{t)),u{t)) 



lim - 

s-»oo S jQ 



'j{u{t),T]u(t){r)),u{t))dT 



< -a\u{t)\'^ -h{riu{t),u{t)) + c 



Then by the Gronwall lemma and ffTUj) there is some positive constant C such 
that for any t > 

E|M(t)p < C{1 + ImoH . 

Now by the same analysis of the proof of Theorem [5] we obtain ( |2^ . Then by 
the above a priori estimates and the local Lipschitz property of /, a standard 
method [12] yields the existence and uniqueness of u. □ 



4 Averaged equation 

This section gives the averaged equation and, as a byproduct, the convergence 
rate is obtained. For this we consider our system in a smaller probability 
space. By Theorem [71 for any k > there is compact set Kf^ in C(0,T; H) 
such that 

Here K^, is chosen as a family of decreasing sets with respect to k. Moreover, 
by the estimate ffTTl) and Markov inequality, we further choose the set K,^ 
such that for u'^ G Kf^ 

\\u^it)f<C^. te[0,T], 
for some positive constant C^. 

Proof of Theorem Now we prove the rate of convergence. In order to do 
this, for any k > we introduce a new sub-probability space {Q^, J-'k,^,^) 
defined by 

Then P(f2 \ Q^) < k. In the following we denote by the expectation 
operator with respect to P^. 

Now we restrict G and introduce an auxiliary process. For any 
T > , partition the interval [0, T] into subintervals of length 5 = i/e . Then 
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we construct processes {u'', v") such that for t G [k6, {k + 1)5), 



kS 



+ (Ti / e^^'-'^dWiis) , u'{0) = uo (25) 



kS 

dv'it) = - [Bd'{t) + g{u'{k6),v'{t))] dt + ^ dW2{t) , 
e ye 

v'(k5) = v'{k5) . (26) 
By the Ito formula for t G [k6, {k + 1)5) 

< -i(Ai - L,)\v-{t} - f.'(i)|' + - £'*(«)l|ii*(i) - "'(M)! . 

By the choice of fi^, is compact in space C{0,T; H), there is > 0, 
such that 

\u'{t) -u'{k5)\^ <C^5^ (27) 

for t G [M, {k + 1)5). Then by the Gronwall lemma, 

\v'{t) -v'{t)\^ <CtS\ te[0,T]. (28) 

Moreover, by the choice of fi^ and the assumption on the growth of f{-,v), 
f{-,v) : H is Lipschitz with —1/2 < /3 < —1/4. Then we have for 

t G [M, (A; + 1)5) 



\u\t) -u'{t)\ < Cf \v'{s)-d'{s)\ds + CfC^ \u'{k6)-u'{s)\ds, 

JkS JkS 

for some positive constant . So by noticing ( 1271) . we have 

\u'{t) -u'{t)\ <C!^6 , tG[0,T]. (29) 

On the other hand, in the mild sense 

t ft 



u{t)=e^'uo+ / e^(*-^)/(M(s))rfs + (Ti / e^^'-'UWi{s) . 
Jo Jo 

Then, using [z\ to denote the largest integer less than or equal to z, 

\u'{t)-u{t)\ < re^(*-)|/(n^(Ls/5J5),F(3))-/V(Ls/5J5))|d3 



*e^(*-)|/V(Ls/5J5))-/V(s))|rfs 







3^(*-)|/V(s)) -/>(.))! 
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Then hj 1^, and we have for t G [0, T] 



\u'{t)-u{t)\ < 



T 



6 + / \u'(s) - u(s)\ ds 







(30) 



As 

\u'{t) - u{t)\ < \u'{t) - u{t) \ + \u{t) - u{t)\ 
by the Gronwall lemma and f l27j) . ( l29l) and ( 130|) we have for t G [0,T], 

|M'(t) -M(t)| < C^rv^. (31) 
The proof of Theorem [2] is complete. □ 

5 Deviation estimate 

The previous section proved that for any T > in the sense of probability 

sup \u'{t) -u{t)\< Cryfe 
*e[o,T] 

for some positive constant Ct- Formally we should have the following form 
M'^(t) = u{t) + (9(e^/^) . This section determines the coefficient of e^^^, the 
deviation. 



Proof of Theorem 0. We approximate the deviation defined by ^ for 
small e > 0. The deviation z"^ satisfies 

z^ = Az^ + ^[f{u^,v^)-f{u)\ (32) 
with z'^{Q) = . By the assumption on / we have 

+ ^{f{u,r^:)-f{u),z^). (33) 
Then the Gronwall lemma yields that for any T > , 

E snp \z'{t)\^ + E [ \\z'{t)fdt<CTi\uo\'^ + \vo\^). (34) 

0<t<T Jo 
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In the mild sense we write 



Then for any < s < t , by the property of e"^*, we have for some positive 
1 > 5 > 



\z'it)-z^is)\ < ^ 



f{u'{r),v'{r)) - f{u{r))dr 



e'"-" ''f{u^{r),v''{r)) — f{u{r))dr 



1 



By the assumption Hi, Theorem [5] and 



E- 



(35) 



Then we have 

nz'it)\cS^O,T;H)<CTi\Uo\'+\Vo\'). 

Here C\0,T;H) is the H51der space with exponent 6. On the other hand, 
also by the property of e"^*, we have for some positive constant CT,a and for 
some 1 > a > 

\z'{t)\H'^ < ^ [\t-s)-"/'\f{u'{s),V%s))-f{u{s))\^ds 

V e Jo 



Then 



E sup \z'it)\H^ <CTAK\'' + \v^o\)■ 
0<t<T 



(36) 



And by the compact embedding of ^-5(0, T; H) f] C(0, T; if") C C(0, T; if), 
{i^''}e , the distribution of {z''}^ , is tight in C(0,T; ii). 
Divide z*^ into zl + which solves 



zl = Azl + -^[f{u,v:) - /»] , ^i(O)^ = 



(37) 



and 



= Az^, + -^[f{u% v^) - f{u, <)] , ^^(0) = 



(38) 



11 



respectively and consider and Z2 separately. We follow a martingale ap- 
proach [71 HH]. Denote by be the probability measure of zl induced on 
space C(0, T; H). For 7 > , denote by UC^H, R) the space of all functions 
from if to M which, together with all Frechet derivatives to order 7, are uni- 
formly continuous. For h G f/C'''(if, M) , denote by h' and h" the first and 
second order Frechet derivative. Then we have the following lemma. 

Lemma 9. Assume H1-H4. Any limiting measure of u^, denote by , 
solves the following martingale problem on C(0,T; H): P^{zi{fS) = 0} = 1 , 



h{zr(t))-h{z^{Q))- I {h\z,{T)),Az,{T))dr-^ I ti [h"{z,{T))B{u)]dT 



is a P° martingale for any h G UC'^{H,'R) . Here 



Biu) 



E[(/(M,77„(t)) - /») ® {f{u,UO)) - /»)] dt 



and ® denotes the tensor product. 

Proof. We follow a martingale approach [TJ [16] . For any < s < t < 00 and 
h G UC^iH) we have 



h{zl{t)) - h{zl{s)) 



{h'{zl{r)),^)dT 



{h'{zl{T)),Azl{T))dT+— I {h'{zl{T))J{u{T),r^l{T))-f{u{T)))dT. 



Rewrite the second term as 



— / {h'{zl{r))J{u{T),v'i{T))-f{u{T)))dT 



— I {h'{zl{t))J{u{T),r^:{T))-f{u{T)))dT 



t ft 



^/ / h"{z{{5)) f{u{r),r^l{T))-f{u{T))),Az\{6) 



S J T 

t rt 



d6 dr 



Li -\- L2 -\- 



d6 dr 
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where Li, L2 and L3 denote the separate hnes of the right-hand side of the 
above equation, respectively. Let {cj}^]^ be one eigenbasis of H, then 

h"{zl{S)) ((/(«(r), r^lir)) - f{u{r))), f{u{S),vm ' 
00 

= J]9.,M^i(5))((/(^(r),r^:(r)) -/Xr))) 

®(/K5),r/:(5))-/K5))),e.®e,>. 

Here dij = de^de where de^ is the directional derivative in direction and $>> 
denotes the tensor product. 

_ Denote by A^,. (5, r) = ((/(«(r), r^^(,)(r)) - /Xr))) ® - 
f {u{S))) , Ci ® Cj) . Then we have 



^3 = --E/ I di,h{z{{6))A\^{6,T)d6dT 

e Js Jr 

-7 E / / / [AzliX) + ^(/(n(A), r^:(,)(A)) - /XA)))] > 



X i-, ((5,r) dXdddr 



+ [ d,M4it))At^iS,T)d6dT 



= L31 + L32 + L33 

with y4-^((5, r) = y4-^ (5, r) — E[74^^.(5, r)]. For our purpose, for any bounded 
continuous function $ on C{0,s;H), let $(-,tu) = $(zj(-,a;)). Then by (dSD, 
we have the following estimate 

|E[(L3i + L32)$]| ^0 ase^O. 

Now we determine the limit of W,AIj{6,t) d6 as e — . Notice that t]u 
depends on u, AIj{6,t) is not a stationary process for fixed r. For this 
introduce 

Z;.(5,r) = ((/(t.(r),r^:(,)(r)) - /Xr))) ® (/(«(r), r/:(,)(5)) - /Xr))), e. ® e,> . 
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Then 



nAU8,T) - AlA5,r)]d5 



< 



E 



-(/Kr),r/:(.)(5))-/Kr)),e,-)] 



d5 



E 



ifHr), <(,)(r)) - /(n(r)), e,) [{f{u{S), ^1^^^ - e,) 



By the assumption Hi and ( 122|) we have 

|(/(«(5),r/:(,)(<5))-/(«(5),r/:(,)(5)),e,)| 



< 



Cf\vtiS)i5)-Vkr)iS)\ 



< CfC\e, 



and by Lemma [H] 



|(/(«W,^:MW)-/(«(r),<(.)(5)),e,)| 



\u{r)r]\u{5) - u{r)\\e,\ 



< Ct{1 



Then by (pS]) 



E[4^,((5,r)-4(<5,r)]rf<5 



, as e — s> . 



Now for fixed r, since riu(r)(t) is stationary correlated, we put 



Then we have 



E[((/(«(r),r/„(.)(5))-/Xr))) 
® (/(^(^),^n(r)(T)) - f{u{r))),ei^ej)] 



r 



Further by the exponential mixing property, for any fixed 6 > t 



{5-r)/e 



e ^ 0. 
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Then, if e„ — > as n — > oo , z/*^" —>■ , 

lim E[L3$] = - / E^°(tT \h"{zi{T))B{u{T))]<^)dT. 

n^oo 2 J g \ / 



where 



Moreover by the assumption on / and the estimates of Lemma [HI B{u) : H — * 
H is Hilbert-Schmidt. 
Similarly by ( |T9l) 

E[Li$ + L2$] ^0 ase^O. 

By the tightness of z" in C{0,T; H), the sequence zl" has a hmit process, 
denote by Zi, in weak sense. Then 



hm E 

n— >oo 



(/i'(z^"(r)),v4z^"(r))$rfr 



E 



{h'{z^{T)),Az,{T))<^dT 



and 



lim E[(/i(2^"(t)) - /i(z^"(s)))$] = E[(/i(zi(t)) - /i(zi(s)))<l>] . 
At last we have 

E'''[{h{z,)it)-h{z^{s)))<t>] (39) 



E 



pO 



pO 



(/i'(^i(r)),A^i(r))<l>c/r 

tr [/;,"(2i(r))5(M(r))]<l>dr 



+ -E 
2 

By an approximation argument we can prove ( l39l) holds for all /i G UC^{H). 
This completes the proof. □ 

We need a lemma on the martingale problem. First introduce some nota- 
tion. Suppose 2t is a generator of bounded analytic compact semigroup S{t). 
d ■■ [0,T] X H ^ H and ^ : [0,T] x H ^ C{H) is B{H) and B{C2{H)) 
measurable and bounded. Let ilt be a second order Kolmogorov diffusion 
operator of the form 

£tF = ^ tr[*Bg^*F,,] + (21^, F,) + F,) 

for any bounded continuous function F on H with first and second order 
Frechet derivatives. Then we have the following result 
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Lemma 10. For any T > , F{X{t)) - F(X(0)) - /J S.sF{X{s)) ds is a P° 
martingale on space C(0, T; H) if and only if the following equation 

dX{t) = [QlX{t) + d{X{t))] dt + ^ dW{t) 

has a weak solution {Q, {^t}o<t<T,^, ^{t),W{t)) such that is the image 
measure o/P hy uo X{-,uj) . 

By the uniqueness of solution, the hmit of z/J, denote by P°, is unique 
and solves the martingale problem related to the following stochastic partial 
differential equation 

dzi = Azi dt + ^B{u) dW , (40) 

where W{t) is cylindrical Wiener process with trace operator Q = Idn, iden- 
tity operator on H, defined on a probability space (0, JF, P) such that zl con- 
verges in probability P to Zi in C(0, T; H). 

On the other hand, the distribution of Z2 on C{0,T;H) is also tight. 
Suppose Z2 is one weak limit point of Z2 in C{0,T; H). We determine the 
equation satisfied by Z2. From (!38|) 



for the convex combination = 9u'' + {l — 9)u , 6* G (0, 1). By assumption Hi 
and ([HD 

-^E\fiu%v^) - /(«^r/:)| < -LcfE\v^-v:\ ^ as e ^ . 
V ^ V ^ 

And for any h G H^, 

E{fi{u%v:X-TMz,h) _ 

— s> as e — s> , 

where z = + 2:2 . Then Z2 solves the following equation 

Z2 = Az2+Ti{u)z, Z2(0)=0. (41) 

And by the wellposedeness of the above problem, we have that z^ uniquely 
converges in distribution to z which solves ( ITUI) . This proves Theorem El □ 
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Figure 1: an example realisation of the FitzHugli-Nagumo system (I42l)-(l43l) 
with cTi = , (72 = 3 , W = (/ — for a cylindrical Wiener process Z(t), 

small parameter e = 0.1 , and domain L = 1 . 

6 Application to stochastic FitzHugh— Nagumo 
system 

Consider the following stochastic FitzHugh-Nagumo system with Dirichlet 
boundary on [—L,L]: 



W{t) is a L'^{—L, L)-valued Wiener process with covariance Q. Let A = d^x 
with zero Dirichlet boundary on {—L,L) , f{u,v) = u — + v , g{u,v) = 

+ M , cTi = and CT2 = 3 , then (jl2D-(jlH!) is in the form of ©-(E]). Figured 
plots an example solution of the FitzHugh-Nagumo system (II1)-(I2]) showing 
that the noise forcing of v feeds indirectly into the dynamics of u. 

Furthermore, for any fixed u the SPDE (1431) has a unique stationary solu- 
tion 1]^ with distribution 



[ul, + u'- {u')^ + v'] dt , 

1 3 

- Hx + u"] dt+— dW. 



(42) 



(43) 



/i„ = AT ( (/ - dxx)-\ Y + - ® " 





(44) 
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Figure 2: RMS of the slow mode u'^{0, t) versus the half-length L at parameter 
e = 0.1 for numerical simulations of (H2l) -( H3l) . This clearly shows the bifur- 
cation that the averaged equation (l44l) predicts from L = 7r/2^/^ = 1.3207 as 
shown by the curve. 
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Figure 3: an example realisation of the stochastic deviation ( H5ll for small 
parameter e = 0.1 , (T2 = 3 , Q = (1 — dxx)~^ and L = 1. 

This averaged equation predicts a bifurcation as L increases. For a fun- 
damental mode on {—L,L) of u = a cos /ex for wavenumber k = 7r/(2L), 
the linear dynamics of the deterministic averaged pde (jH]) predicts that 
M = is stable for k > 2^/'^, that is, L < For larger domains with 

L > 7r/2^/^ = 1.3209 the averaged pde predicts a bifurcation to finite am- 
plitude solutions. This bifurcation matches well with numerical solutions as 
seen in Figure [2] which plots the mean mid-value as a function of L: the 
bifurcation is clear albeit stochastic. 

To quantify the fluctuations evident in the dynamics of u'^ we turn to the 
PDE for deviations. Denote by rju the unique stationary solution of 

dv = [dxxV — V + u] dt + dW. 

Then f{u, rju) — f{u) = rjuit) — (/ — dxx)~^u and we have that the deviation z 
solves the spde 

dz = [dxxz + (1 - Su'^)z] dt + 3(/ - d^x^^^/QdW (45) 

with W{t) being a cylindrical Wiener process defined on a larger probability 
space with covariance operator Id on L^([— L, L]). Figure [3] plots an example 
of the deviation between the original system and the averaged system, the 
SPDE ( l45l) . Including the deviation spde ( l45l) gives a much better approxi- 
mation than the deterministic averaged equation (jBI). In particular, when 
the initial state Uq = and there is no direct forcing of m, ai = , as used in 
Figures [T] and El then the averaged solution is identically u{t) = . In such a 
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Figure 4: simulations over a time of 128 (with L = 1) show: circles, fluctu- 
ations in u^{0,t) whose variance is plotted as function of e {At — 0.002); 
crosses, fluctuations in ^/ez whose variance is plotted as a function of e 
(At = 0.02). The two agree remarkably well: their standard deviations 
are best fitted by the nearly identical O.OGSy^ and 0.066v^, respectively. 
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case, the dynamics of u as seen in Figure [T] are modelled solely by deviations 
governed by the SPDE 

To quantify the comparison between the deviation pde f HSl) and the orig- 
inal dynamics of the FitzHugh-Nagumo system (l42l)-(H3l). we look at how 
the fluctuations scale with scale separation parameter e. At L = 1 the noise 
free state = = is stable. As parameter e increases the fluctuations 
in have variance as plotted by circles in Figure HI The scatter in the plot 
reflects that averages over much longer times would be better. However, the 
variance does scale with e as required, and in close correspondence to that 
predicted by the deviation pde ( H5l) (crosses). 
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